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Abstract 

The density matrix of composite spin system is discussed in relation to the 
adjoint representation of unitary group U(n). The entanglement structure 
is introduced as an additional ingredient to the description of the linear 
space carrying the adjoint representation. Positive maps of density operator 
are related to random matrices. The tomographic probability description of 
quantum states is used to formulate the problem of separability and entan- 
glement as the condition for joint probability distribution of several random 
variables represented as the convex sum of products of probabilities of ran- 
dom variables describing the subsystems. The property is discussed as a 
possible criterion for separability or entanglement. The convenient criterion 
of positivity of finite and infinite matrix is obtained. The [/(n)-tomogram 
of a multiparticle spin state is introduced. The entanglement measure is 
considered in terms of this tomogram. 

KEY WORDS: unitary group, entanglement, adjoint representation, to- 
mogram, operator symbol, random matrix. 



1 Introductiion 

The notion of entanglement pQ is related to the quantum composition prin- 
ciple of the states of subsystems for a given multipartite system. For pure 
states, the notion of entanglement and separability can be given as follows. 

If the wave function of a state of a bipartite system is represented as the 
product of two wave functions depending on coordinates of the subsystems, 
the state is simply separable; correspondingly, in other cases, the state is en- 
tangled. An intrinsic approach to the entanglement measure was suggested 



1 



in j2] . The measure was introduced as the distance between the system den- 
sity matrix and the tensor product of the subsystem states. There are several 
other different characteristics and measures of entanglement considered by 
several authors [3-9]. Each of the entanglement measures describes a degree 
of correlations between the subsystems' properties. The notion of entan- 
glement is not an absolute notion for a given system but depends on the 
decomposition into subsystems. The same quantum state can be considered 
as entangled, if one kind of division of the system into subsystems is given, 
or as completely disentangled, if another decomposition of the system into 
subsystems is considered. 

For instance, the state of two continuous quadratures can be entangled 
in Cartesian coordinates and disentangled in polar coordinates. Coordinates 
are considered as measurable observables labeling the subsystems of the given 
system. The choice of different subsystems mathematically implies the exis- 
tence of two different sets of the subsystems' characteristics (we focus on bi- 
partite case). We may consider the Hilbert space of states H(l, 2) or H(l', 2'). 
The Hilbert space for the total system is, of course, the same but the index 
(1,2) means that there are two sets of operators P\ and P 2 , which select 
subsystem states 1 and 2. The index (l',2') means that there are other two 
sets of operators P[ and P 2 , which select subsystem states 1' and 2'. The 
operators P\p and P[, 2 , have specific properties. They are represented as 
tensor products of operators acting in the space of states of the subsystem 1 
(or 2) and unit operators acting in the subsystem 2 (or 1). In other words, 
we consider the space H, which can be treated as tensor product of spaces 
H(l) and H(2) or H(V) and H(2'). In the subsystems 1 and 2, there are 
basis vectors | ni) and | m 2 ), as well as in the subsystems 1' and 2' there 
are basis vectors | n[) and | m' 2 ). The vectors | ri\) \ m 2 ) and the vectors 
I n '\) I m 2) f° rm the sets of basis vectors in the composite Hilbert space, 
respectively. These two sets are related by means of unitary transformation. 
An example of such a composite system is a bipartite spin system. 

If one has spin-ji [the space H(l)] and spin-j 2 [the space H(2)\ systems, 
the combined system can be treated as having basis | jimi) \ j 2 m 2 ). 

Another basis in the composite-system-state space can be considered in 
the form | jm), where j is one of the numbers \j\ — j 2 \, \j\ — j 2 \ + 1, . . . , 
ji + J2 and m = mi + m 2 . The basis | jm) is related to the basis | j±mi) | 
j 2 m 2 ) by means of unitary transform given by Clebsch-Gordon coefficients 
C {j\Wi\j 2 m 2 \jm) . From the viewpoint of given definition, the states | jm) 
are entangled states. For example, if j\ = j 2 = 1/2, there are entangled spin 
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states of the composite system, which nowadays are called Bell states 



These states are maximally entangled states. In terms of spin, the states 
| $ ± ) are the superpositions of j — 1, m — ±1 states and the states | 1 I r± ) 
are the superpositions of j = 0, 1, m — states. 

The spin states can be described by means of the tomographic map [10- 
12]. For bipartite spin systems, the states were described by the tomographic 
probabilities in f3J • Some properties of the tomographic spin description 
were studied in [15| . In the tomographic approach, the problems of the quan- 
tum state entanglement can be cast into the form of some relations among 
the probability distribution functions. On the other hand, to have a clear 
picture of entanglement, one needs mathematical formulation of properties of 
the density matrix of the composite system, a description of the linear space 
of the composite system states. Since the density matrix is hermitian, the 
space of states is a subset of linear space of adjoint representation of the group 
U(n 2 ), where n = 2j + 1 is the dimension of the spin states of two spinning 
particles. Thus one needs to characterize the connection of the entanglement 
phenomena with the structures in the space of adjoint representation of the 
U{n 2 ) group. 

The aim of this paper is to connect entanglement problems with the 
properties of tomographic probability distributions and discuss the properties 
of the convex set of positive states for composite system with taking into 
account the subsystem structures. We used Hilbert-Schmidt distance to 
calculate the measure of entanglement as the distance between a given state 
and the tensor product of the partial traces of the density matrix of the 
given state. In jTHj another measure of entanglement as a characteristic 
of subsystem correlations was introduced. This measure is determined via 
covariance matrix of some observables. Review of different approaches to the 
entanglement notion and entanglement measures is given in jT7j, where the 
approach to describe entanglement and separability of composite systems is 
based on entropy methods. 

Due to variety of approaches to the entanglement problem, one needs 
to understand better what in reality this word 'entanglement' describes. Is 





3 



it a synonym of the word 'correlation' between two subsystems or it has 
to capture some specific correlations attributed completely and only to the 
quantum domain? 

The paper is organized as follows. 

In Sec. 2 we study the subsystem structure of a given linear space. In Sec. 
3 we consider the relation of the group U (n 2 ) to the set of density matrices. In 
Sec. 4 we discuss positive maps. In Sec. 5 we investigate local transforms. In 
Sec. 6 we treat the probability distributions as vectors. In. Sec. 7 we prove 
the invariance of the intrinsic entanglement measure. In Sec. 8 we define the 
separable states. In Sec. 9 generic symbols of operators are presented. In 
Sec. 10 an example of Weyl symbols is considered and in Sec. 11 an example 
of quadrature tomogram is done. In Sec. 12 symbols of density operators for 
multipartite system are discussed. Spin tomography is reviewed in Sec. 13. 
Two qubits are considered in the tomographic representation in Sec. 14. The 
relation of dynamical map to purification procedure is described in Sec. 15. 
Some properties of quadratic forms are reviewed in Sec. 16. The tomogram 
for the group U(N) is introduced in Sec. 17. Conclusions and results are 
listed in Sec. 18. 

2 Linear Space of a Composite System, Its 
Structure, and Its Convex Subset of Posi- 
tive States 

In this section, we review the meaning and notion of composite system in 
terms of additional structures on the linear space of state for the composite 
system. 

2.1 States and Observables 

In quantum mechanics, there are two principal ingredients, which are asso- 
ciated with linear operators acting in a Hilbert space. The first ingredient is 
related to the concept of quantum state and the second one, to the concept of 
observable. The state is associated to Hermitian nonnegative, trace-class, lin- 
ear operator. The observables are associated to Hermitian operators. Though 
the both states and observables are identified with the Hermitian operators, 
there is an essential difference between these two objects. The observables 
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have additional product structure. Thus we consider product of two linear 
Hermitian operators corresponding to the observables. First measuring an 
observable and (after measuring the first one) measuring another observable 
just correspond to the product of two operators. 

For the states, the notion of product is redundant. The product of two 
states is not a state. For states, one keeps only the linear structure of vector 
space. For finite n-dimensional system, the Hermitian states and the Her- 
mitian observables live in Lie algebra of the unitary group U(n). But the 
states correspond to nonnegative Hermitian operator. The observables can 
be associated with both types of the operators including nonnegative and 
nonpositive ones. Space of states is linear space which, in principle, is not 
equipped by a product structure. Due to this, if one considers transforma- 
tions in linear space of states, one does not need to preserve any product 
structure. In the set of observables, one needs to care what is happening 
with product of operators provided some transformations are applied. 

2.2 Vectors 

Let us first introduce some extra constructions of the map of a matrix onto 
a vector. Given a rectangular matrix M with elements M^, where i = 
1,2, ... ,n and d — 1,2, ... ,m. Then one can consider the matrix as a vector 
Ai with N = nm components constructed by the following rule: 

Mt = Mn, M 2 = M 12 , M m = M lm , M m+l = M 21 , . . . M N = M nm . 

Thus we construct the map M — > M. =tj^ M M. 

We have introduced the linear operator t^ M which maps the matrix 
M on a vector Ai. Now we introduce the inverse operator pj^ M which 
maps a given vector column in the space with dimension N = mn onto the 
rectangular matrix. This means that given a vector M. = Ai±, . . . , A4n, we 
use the rule of relabeling the components of the vector introducing two indices 
i = 1, . . . , n and d = 1, . . . , m. The relabeling is accomplished according to 
(PJ. Then we collect the relabeled components into matrix table. Thus we 
get the map 

p MM M = M. (2) 
One can see that their composition 

ImmPAm-M =1 ■ M (3) 
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acts on the vector as unit operator in the linear space of vectors. 

Given anxii matrix the map suggested can also be extended. The ma- 
trix can be treated as n 2 -dimensional vector and, vice versa, the vector of 
dimension n 2 can be mapped by this procedure onto the nxn matrix. 

— * 

Let us consider a linear operator acting on the vector Ai and related to 
a linear transform of the matrix M. First, we study the correspondence of 
the linear transform of the form 

M -> gM = M l g (4) 

to the transform of the vector 

M -> M\ = C l g M. (5) 

One can show that the n 2 xn 2 matrix C l g is determined by the tensor product 
of the nxn matrix g and nxn unit matrix, i.e., 

4 = <7®1. (6) 

Analogously, the linear transform of the matrix M of the form 

M — > Mg = M r g (7) 

— * 

induces the linear transform of the vector M. of the form 

M^M r g = i MM M r g =ClM, (8) 
where the n 2 xn 2 matrix C r g reads 

C g = l® g tT . (9) 
Similarity transformation of the matrix M of the form 

M - gMg~ l (10) 
induces the corresponding linear transform of the vector A4 of the form 

M^M, = £ 9 g M, (11) 
where the n 2 xn 2 matrix C g reads 

-i\tr 



C s q = g®(g-T- (12) 
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One can ask how to determine the inverse map of vector M onto matrix 
M, i.e., how to define the operator Vmm- ^ n f &c ^ the reconstruction can be 
defined by means of star-product of vectors All in a linear space. One can 
define the associative product of two N- vectors Mi and M2 using the rule 

M = M 1 *M 2 , (13) 

where 

M k = £ Kt(MMM 2 ) s . (14) 

l,s=l 

If one applies a linear transform to the vectors Mi, Mi, M of the form 

Mi -> aI; = A4 2 -> A4' 2 = ^Al 2 , A* -> A? = £A4, 

the invariance of the star-product kernel yields 

Af'^Alf, = Af', if C = G®G- ltT , GeGL{n). 

The kernel If^ (structure constants) which determines the associative star- 
product satisfies the quadratic equation. Thus if one wants to make the cor- 
respondence of the vector star-product to the standard matrix product (row 
by column), the matrix M must be constructed appropriately. For example, 
if the vector star-product is commutative, the matrix M corresponding to 
the iV-vector M can be chosen as diagonal NxN matrix. This consider- 
ation shows that the map of matrices on the vectors provides star-product 
of the vectors (defines the structure constants or the kernel of star-product) 
and, conversely, if one has the vectors, the map of the vectors onto the ma- 
trices with the standard multiplication rule is determined by the structure 
constants (or by the kernel of the vector star-product). 

The constructed map of matrices on the vectors gives a possibility to en- 
large the dimensionality of the group acting in the linear space of matrices in 
comparison with the standard one. Thus, given a nxn matrix M the left ac- 
tion, the right action, and similarity transformation of the matrix are related 
to the complex group GL(n). On the other hand, the linear transformations 
in the linear space of n 2 -vectors M obtained by using the introduced map 
are determined by the matrices belonging to the group GL(n 2 ). There are 
transformations on the vectors which cannot be simply represented on ma- 
trices. If M — > $(M) is a linear homogeneous function of the matrix M, we 
may represent it by 

&ab = B aa ^ W M a i V . 
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Under rather clear conditions, B aa i bb i can be expressed in terms of its non- 
normalized left and right eigenvectors: 

B aa >,bb> = Y, X aa>(v)ylb>( U )i 
v 

which corresponds to 

n 2 

$(M) = xiWV = J2 x(v)My\v). 

There are possible linear transforms on the matrices and corresponding 
linear transforms on the induced vector space which belong not to a group 
but to an algebra of matrices. One can describe the map of nxn matrices M 
(source space) onto vectors M. (target space) using specific basis in the space 
of the matrices. The basis is given by the matrices Ej k (j,k = 1,2, ... ,n) 
with all matrix elements equal to zero except the element in jth row and kth 
column which is equal to unity. One has the obvious property 

M jk = Tr(ME jk ). (15) 

In our procedure, the basis matrix Ejk is mapped onto the basis column- 

— * 

vector Ejk, which has all components equal to zero except the unity compo- 
nent related to the position in the matrix determined by the numbers j and 
k. Then one has 

n 

M= Tr(ME jk )S jk . (16) 

j,k=i 

For example, for similarity transformation of the finite matrix M, one has 

N 

M" g = E Tr (gMg-'E^) £ jk . (17) 
j,k=i 

Now we will define the notion of 'composite' vector which corresponds to 
dividing a quantum system into subsystems. 
We will use the following terminology. 

In general, the given linear space of dimensionality N = van has a struc- 
ture of bipartite system, if the space is equipped with the operator pj^ M 
and the matrix M (obtained by means of the map) has matrix elements in 
factorizable form 

M id -> x t y d . (18) 



8 



This M = x <g> y corresponds to the special case of nonentangled states. 
Otherwise, one needs 

M = J2x(u)(g)y(iy). 

V 

In fact, to consider in detail the entanglement phenomenon, in the bipartite 
system of spin, one has to introduce a hierarchy of three linear spaces. The 
first space of pure spin states is two-dimensional linear space of complex 
vectors 

(:;)• 



X ) 



(19) 



In this space, the scalar product is defined as follows: 

(x\y) = x\yi +x* 2 y 2 . 



(20) 



So it is two-dimensional Hilbert space. We do not equip this space with a 
vector star-product structure. In the primary linear space, one introduces 
linear operators M which are described by 2x2 matrices M. Due to the map 
discussed in the previous section, the matrices are represented by 4-vectors 
M belonging to the second complex 4-dimensional space. Star-product of 
the vectors M determined by the kernel JCf s is defined in such a manner in 
order to correspond to the standard rule of multiplication of the matrices. 

In addition to the star-product structure, we introduce the scalar product 
of the vectors Mi and M 2 , in view of the definition 



(Mi | M 2 ) =Tr (MlM 2 ), 



(21) 



which is the trace formula for scalar product of matrices. 

This means introducing the metric g a/3 in the standard notation for scalar 
product 



(All I M 2 ) = £ (MO^Ma),, 



(22) 



where the matrix g alS is of the form 



( 1 








\ 








1 








1 








V o 








1 J 



9 J 9 



(23) 
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The scalar product is invariant under action of the group of nonsingular 4x4 
matrices £, which satisfy the condition 

r 1 = gtg. (24) 

The product of matrices I satisfies the same condition since g 2 = 1. 

Thus, the space of operators M in primary two-dimensional space of 
spin states is mapped onto linear space which is equipped with a scalar 
product (metric Hilbert space structure) and associative star-product (kernel 
satisfying quadratic associativity equation). In the linear space of the 4- 

— # 

vectors Ai, we introduce linear operators (superoperators), which can be 
associated with algebra of 4x4 complex matrices. 

3 Density Operators and Positive Maps 

In this section, we focus on density matrices. This means that our matrix 
M is considered as a density matrix p which describes a quantum state. We 
consider here the action of the unitary transformation U(n) of the density 
matrices and corresponding transformations on the vector space. If one has 
structure of bipartite system, we also consider the action of local gauge trans- 
formation both in the 'source space' of density matrices and in the 'target 
space' of the corresponding vectors. 

The nxn density matrix p has matrix elements 

Pik = pl, Trp = l, <##)>0. (25) 

Since the density matrix is hermitian, it can be always identified as an element 
of the convex subset of the linear space associated with the Lie algebra of 
U (n) group, on which the group U (n) acts with the adjoint representation 

p^ Pu = UpUl (26) 

The system is said to be bipartite, if the space of representation is equipped 
with an additional structure. It means that for 

n 2 = n 1 - n 2 , n\ = n 2 = n 

one can make first the map of nxn matrix p onto n 2 -dimensional vector p 
according to the previous procedure, i.e., one equips the space by an operator 



10 



tpp. Given this vector one makes a relabeling of the vector p components 
according to the rule 

P^Pid,ke> = 1,2, ...,m, d,e= 1,2, ...,71-2, (27) 
i.e., obtaining the quadratic matrix 

Pq = V Pq pP- (28) 

The unitary transform (j2Ej) of the density matrix induces the linear transform 
of the vector p of the form 

p^Pu = {U®U*)p. (29) 

There exist linear transforms (called positive maps) of the density matrix, 
which preserve its trace, hermicity, and positivity. It is the transform intro- 
duced in [18j 

Po ^ p s = J2PkUkpoUt, X)P* = 1 > (30) 

k k 

where Uk are unitary matrices and pk are positive numbers. 

If the initial density matrix is diagonal, i.e., it belongs to Cartan sub- 
algebra of Lie algebra of the unitary group, the diagonal elements of the 
obtained matrix give smoother probability distribution than the initial one. 
There exists the generic transform (see )T%) 119j) 

Po - P = EVWtf, Y,VkVk = l- (31) 

k k 

For large number of terms in the sum, the above map gives the most stochas- 
tic density matrix 

Po -> Ps = (ny l \. 

The transform (j30j) is the partial case of the transform 1)31)1 . We discuss the 
transforms separately since they are used in the literature in the presented 
form. 

One can see that the constructed map of density matrices onto vectors 
provides the corresponding transforms of the vectors, i.e., 

Po -> Ps =J2Pk( U k® U k)Po (32) 

k 
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and 

Po^P = £(^®K)A). (33) 

k 

It is obvious that the linear transforms of the vectors, which preserve their 
properties to correspond to the density matrix, are essentially larger than 
the standard unitary transform of the density matrix. 

Formulas and mean that the positive map superoperators acting 
on the density matrix in the vector representation are described by n 2 xn 2 
matrices 

C s = J2Pk(Uk ® Ul) (34) 

k 

and 

£ = J2V k ®V k *, (35) 

k 

respectively. 

Positive map is called 'noncompletely positive' if 

c = e v h ® v k * - y: v s ® e vtfn - e «t w- = i. 

fe s A; s 

This map is related to nonphysical evolution of a subsystem. 



4 Positive Map and Random Matrices 

Formula ()34j) can be considered in the context of random matrix representa- 
tion. In fact, the matrix C s can be interpreted as the weighted mean value 
of the random matrix Uk®U£. The dependence of matrix elements and pos- 
itive numbers pk on index k means that we have a probability distribution 
function pk and averaging of the random matrix Uk <8> U% by means of the 
distribution function. So the matrix C s reads 

C S = (U®U*). (36) 

Let us consider example of 2x2 unitary matrix. We can consider the matrix 
of SU (2) group of the form 

A \a\ 2 +\(3\ 2 = l. (37) 
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The 4x4 matrix C s takes the form 



—n 
—n* 
\l-t 



m 
s 



-q 

-m 



m 

-Q 

s* 
-m* 



n 
n* 



The matrix elements of the matrix C s are the means 

m = (a(3*), 

£ = (aa*), 

n = (aft), 

s = (a 2 ), 

Q = (P 2 )- 

Moduli of these matrix elements are smaller than unity. 
Determinant of the matrix £„ reads 



det C s 



2*) (I 



4 Re 



q m n + vans 



If one represents the matrix C s in block form 



then 



and 



A 



t m 
—n s 



A B 
C D 



B 



m 



D = a 2 A*a 2 , 



V —q n 

C = —a 2 B*a 2 , 



(38) 



(39) 



(40) 



(41) 



(42) 



(43) 



where a 2 is Pauli matrix. 

One can check that the product of two different matrices C s can be cast 
in the same form. This means that the matrices C s form the 9-parameter 
compact semigroup. For example, in the case i = 1/2 and m — 0, one has 
the matrices 

1/2 ^ B = ( 1/2 ^ 



A = 



-n s 



-q n 



(44) 
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Determinant of the matrix C s in this case is equal to zero. All the matrices 
C s have the eigenvector 

( l l 2 \ 




V 1/2 / 



Po 



(45) 



i.e. 



This eigenvector corresponds to the density matrix 



Pi 



1/2 
1/2 



(46) 



(47) 



which is obviously invariant of the positive map. 

For random matrix, one has correlations of the random matrix elements, 
e.g., (aa*) ^ (a)(a*). 



The matrix C r 



maps the vector 



^ 1 ^ 

10 

10 

V 1 / 



fPn \ 

Pl2 
P21 
V P22 J 



(48) 



(49) 



onto the vector 



Pt 



(50) 



( Pn \ 

P21 
Pl2 
V P22 J 

This means that the positive map (|4*K|) connects the positive density matrix 
with its transposed. This map can be presented as the connection of the 
matrix p with its transposed of the form 



P^ P 



T 



P 



-z(p + oipoi - V2PV2 + 03/00-3) . 
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There is no unitary transform connecting these matrices. 

This noncompletely positive map in iV-dimensional case is given by gen- 
eralized formula 

1 + e 

p—>p s = -ep + -j-j- 1 N , e > 0. 

The standard unitary transform can be interpreted as average random trans- 
form with probability distribution 

Pk = 6(k), (51) 

where 8{k) is either Kronecker symbol Sko for discrete index k or Dirac delta- 
function for continuous index k. 

For standard unitary transform, one cannot find the matrix U satisfying 
the equation 

U®U* = C p . (52) 

But if one makes averaging with generic distribution function [not with prob- 
ability distribution (pTTj)]. the equation 

(U®U*) = C P (53) 

has the solution. 

The standard unitary transform of density matrix is 3-parameter subset 
of this 9-parameter semigroup. 

Thus we constructed matrix representation of positive map of density 
operators of spin- 1/2 system. To construct this representation, one needs 
to use the map of matrices on the vectors discussed in the previous section. 
Formulas (|31|) and (J35|) can be interpreted also in the context of random 
matrix representation, but we use the uniform distribution for averaging in 
this case. So one has equality (|33j) in the form 

C = (V V*) (54) 

and the equality 

(VV) = 1, (55) 

which provides constrains for used random matrices V. 

Using random matrix formalism, the positive (but not completely posi- 
tive) maps can be presented in the form 

C=(V®V*)-(v®v*), (V ] V) - (v ] v) = 1. 
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In [IB] the positive maps ([50)1 and (|3T|) were used to describe non-Hamiltonian 
evolution of quantum states for open systems. If the map corresponds to 
an extended Hamiltonian evolution, the leading terms are of order t 2 (Zeno 
effect) and consequently, in the A 2 t-approximation of Prigogine and van Hove, 
one can derive rate equations from macroscopic equations of motion. 

We have to point out that, in general, such evolution is not described by 
first-order-in-time partial differential equation. Like in the previous case, if 
there are added structures of the matrix in the form 



Pid,ke -> Xiy d Z k t e , (56) 

which means association with the initial linear space two extra linear spaces in 
which Xi, Zk are considered as vector components in the ni-dimensional linear 
space and yd, t e are vector components in ^-dimensional vector space, we 
will tell that one has bipartite structure of the initial space of state [bipartite 
structure of the space of adjoint representation of the group U{n)]. Usually 
the adjoint representation of any group is defined per se without any reference 
to possible substructures. Here we introduce the space with extra structure. 
In addition to be space of adjoint representation of the group U(n), it has 
structure of bipartite system. The generalization to multipartite (iV-partite) 
structure is straightforward. One needs only the representation of positive 
integer n 2 in the form 

N 

n 2 = J] nl (57) 
k=i 

If one considers more general map given by superoperator (|35|1 rewritten 
in the form 

c = (v®v*), (y V) = i, 

the number of parameters determining the matrix C can be easily evaluated. 
For example, if 



V = , , V* 





where matrix elements are the complex numbers, the normalization condition 
provides 4 constrains for the real and imaginary parts of matrix elements of 
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the following matrix: 



C 



I (\a\ 2 ) (ab*) (ba*) (bb*) \ 

(ac*) (aa*) (be*) (bd*) 

(ca*) (cb*) (da*) (db*) 

V (cc*) (cd*) (dc*) (dd*) J 



namely, 



<|a| 2 > + (|&| 2 > = l, 



(\c\ 2 ) + (\d\ 2 ) = 1, 



(a*c) + (a*d) = 0. 



Due to structure of the matrix C, there are 6 complex parameters 



{ab*), (ac*), (ad*), (be*), (bd*), (cd*) 



or 12 real parameters. 

Geometrical picture of positive map can be clarified if one considers trans- 
form of the positive density matrix onto another density matrix as transform 
of ellipsoid into another ellipsoid. The generic positive transform means a 
generic transform of the ellipsoid, which changes its orientation, values of 
semiaxis, and position in the space. But the transform is not making from 
the ellipsoid the surface like hyperboloid or paraboloid. For pure states, the 
positive density matrix defines the quadratic form which is maximally de- 
generated. In this sense, we say "ellipsoid" also including all its degenerate 
forms corresponding to density matrix of ranks less than n (in n-dimensional 
case). The number of parameters defining the map (V<g>V*) in n-dimensional 
case is equal to n 2 (n 2 — 1). 



In this section, we discuss the transforms of the density operators of bipartite 
system. We concentrate on the case of two spin- 1/2 systems. Let the first 
system be in the state with 2x2 density matrix 



5 Local and Nonlocal Transforms 




(58) 



and the second system is in the state with 2x2 density matrix 




(59) 
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According to the suggested map, we associate with the state density matrices 
(|58|) and (J59|) the 4-vectors with components 



PA 



fpfi 

Pu 

P21 
V P22 



PB 



{P?l 

Pfl 
V P22 



(60) 



The density matrices (J58j) and (JoTJj) [i.e., vectors ([SUjl] belong to linear spaces 
of adjoint representations of groups 6^4(2) and Ub{2), respectively. 

For the product state (simply separable state) of composite system with 
4x4 density matrix 

Pab = Pa®Pb, (61) 
the corresponding 16- vector associated to 4x4 density matrix 



Pab 



has the form 



( PnPn PnP?2 Pf 2 Pn PX2PX2 \ 

P11P21 P11P22 P12P21 P12P22 

P21P11 P21P12 P22P11 P22P12 

V P21P21 P21P22 P22P21 P22P22 / 



Pab 



(62) 



Pb 



Pa®Pb = Cp An- 



nexe the 16x16 matrix C acting on the 16-component vector p^i 
is standard tensor-product of two vectors, has the form 



(63) 
1 p*B , which 



c 



( h 




















\ 








1 2 




















h 





























1 2 


























1 2 





























1 2 




















1 2 








I 




















h J 



(64) 



The matrix C consists of 2x2-block zero and unity matrices. 

The linear space of Hermitian matrices is equipped also by commutator 
structure defining Lie algebra of the group U(n). The kernel, which defines 
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this structure (Lie product structure) is determined by the kernel, which 
determines star-product. 

In the space of 16-vectors, one defines the scalar product as follows: 

16 

P\-P2= Pla C »/3p2f3- (65) 

a,/3=l 

The product is invariant, in view of linear transform (group transform) 

Pi -> A = Lpi p2 -> P 2 = L P2, (66) 
if the 16x16 matrix L satisfies the condition 

L- 1 = CtfC. (67) 

Vector (JfiHj) belongs to linear space of adjoint representation of the unitary- 
group U (4) and star-product of the vectors is identified with the associative 
algebra generated by the Lie algebra. The local gauge transformations are 
defined as tensor product of independent unitary transforms of the density 
matrices pa ((SHj) an d pn dSSI - These transformations are described by the 
group £7(2) x £7(2). In the space of vectors pA and ps, the local transform 
superoperators have the matrix form 

U A {2)pA = (U A {2) <g> U* A (2))p A (68) 

and 

Ub{2)pb = (U B (2) ® U* B {2))p B . (69) 

In the 16- dimensional space of vectors p A B, the local transforms are described 
by the superoperator 

PAB^ti C B=£ l ° C PAB, (70) 

with the matrix 

C loc = C(U A {2) ® £^(2)) ® (U B (2) ® U B (2))C. (71) 

The local positive map transforms induce in the space of adjoint represen- 
tation of the group U (4) the transform of the vectors pab associated to the 
matrices 

£pu = CfepuUf ® Ut) ® (E^'^ ® U k*)C (72) 
k k' 
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and 

£ P v = C(Y,V k A ® V k M ) ® (£V£ ® V$*)C, (73) 
k k' 

respectively. 

The matrix (J72|) can be expressed in terms of semigroup matrices La and 
Lb as follows: 

L pU = C(L A ®L B )C, (74) 

where 

L A ={U A ®U A *), L B = {U B ®U B *). (75) 

Analogously 

L pV = c({V A ®V M )®(V B ®V B *))C. (76) 
The matrices L v y form 18-parameter semigroup. 

6 Distributions as Vectors 

The notion of entanglement can be better clarified using the concept of dis- 
tance between the quantum states. In this section, we consider the notion of 
distance between the quantum states in terms of vectors. First, let us discuss 
the notion of distance between conventional probability distributions. This 
notion is well known in the classical probability theory. 

Given probability distribution P(k), k = 1,2,...N, one can introduce 
vector P in the form of column with components Pi = -P(l), Pi = -P(2), . . . , 
P N = P(N). The vector satisfies the condition 

N 

E^ = l- (77) 

k=l 

The set of the vectors does not form a linear space but only a convex subset. 
Nevertheless, in this set one can introduce distance between two distributions 
using the vector intuition 

D 2 = (P x - P 2 ) 2 = E p ikPik + E p 2kP2k - 2 E p ikP%k. (78) 

k k k 

One can use another identification of distribution with vectors. 
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Since all P(k) > 0, one can introduce Vk = \JP(k) as components of 

vector V. The V can be thought as column with nonnegative components. 
Then the distance between the two distributions takes the form 

V 2 = (ft - V 2 ) 2 = 2 -2^ v /P 1 (fc)P 2 (fc). (79) 

k 

Two different definitions ()77|) and ()78|) can be used for the notion of 
distance between the distributions. 

Let us discuss now the notion of distance between the quantum states 
determined by density matrices. In the density-matrix space (in the set 
of linear space of adjoint U(n) representation), one can introduce distances 
analogously. The first case is 

Tr( Pl -p 2 ) 2 = D 2 (80) 

and the second case is 

Tr(VpT- Jp- 2 f = V 2 . (81) 

In fact, the distances introduced can be written naturally as norms of vectors 
associated to density matrices 

D 2 = | Pl 

and 

v 2 = ((Vp,) 

respectively. 

In the above expressions, we use scalar product of vectors p\ and p 2 as 
well as scalar products of vectors (VPi) and (y/p 2 ), respectively. 

Both definitions immediately follow by identification of matrices either 
pi and p 2 with vectors according to the map of the previous sections or 
matrices V~P~i and VP~2 with vectors. Since the density matrices p\ and p 2 
have nonnegative eigenvalues, the matrices Jp~{ and VP~2 are defined without 

ambiguity. This means that the vectors (y/p^ and (y/p 2 ) are also defined 
without ambiguity. 

One can easily see that the distance given as the norm of real vector is in- 
variant of the orthogonal group including the improper transform, which are 
discrete transforms like permutations of the vector components and changing 
sign of the components. The norm of the vector is invariant with respect to 
all these transforms. Thus, the invariance group of the norm contains the 
standard local rotations of the orthogonal group and the discrete transforms. 



P2\ 



-(VP: 



2> ' 



(82) 
(83) 
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7 Invariance of Intrinsic Measure 

Given the density matrix pab in the linear space equipped with bipartite 
structure. Then the partial traces exist 

p A = TtbPab, pb = Tt a Pab- (84) 

The Hilbert-Schmidt distance between pab and pa ® Pb, i-e., 

e = Tr (pab - Pa® Pb) 2 (85) 

was considered to define the parameter e as intrinsic measure of the state 
entanglement. It was not proved that this measure is invariant under local 
group U(n\)xU(n<2) transformations. The proof is straightforward. One uses 
the structure of the tensor product in the form 



U nin2 = U(m) ®U(n 2 ) 



( U U V U 12 V . . . U ln V \ 

u 21 v u 22 v ... u 2n v 



(86) 



V UmV u n2 v . . . u nn v j 

where Un~ are matrix elements of the U(n{) group and V is the matrix of 
U{n 2 ) group. 

Due to unitarity 

tC 2 = U j nin2 (87) 
and using the form of this matrix given by Eq. (|86|). one has 



Tr B (U nin2 p AB Ul in2 ) = p A = Uin^pAU^inx), (88) 

Tr^ (U nin2 pA B Ul n2 ) =p B = Vp B Vl (89) 

This means that 

Tr (pab - Pa® Pb) 2 = Tr (p AB ~ Pa® Pb) 2 ■ (90) 

Thus the entanglement e is invariant under local transformations. One can 
introduce another measure of entanglement, which is also invariant under 
local transformations 

1/2 , 1/2 , 1/2 2 

e - Pab - tr A Pab tT B Pab ■ 

Since we have shown that norm of a vector in adjoint representation of the 
unitary group is invariant under action of the group, which is larger than the 
local group, the measures of entanglement introduced are invariant under 
action of local transformations. 
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8 Separable Systems and Separability Crite- 
rion 



According to known definition, the system density matrix is called separable 
(for composite system) if one has decomposition of the form 

PAB = Y.Pk{pA®PB ) ), E^ = 1 > l>Pk>0. (91) 

k k 

The formula does not demand orthogonality of the density operators and 
p^ for different k. Since every density matrix is a convex set of pure density 
matrices, one could demand that pjp and p^ be pure. This formula can be 
interpreted in the context of random matrix representation. In fact, one has 

Pab = (pA® Pb), (92) 

where pa and pB are considered as random density matrices of the subsystems 
A and B, respectively. 

There are several criteria for the system to be separable. We suggest in the 
next sections a new approach to the problem of separability and entanglement 
based on the tomographic probability description of quantum states. The 
states which cannot be represented in the form (|91|) by definition are called 
entangled states [T7j. Thus the states are entangled if in formula (}9~T]) at least 
one coefficient (or more) Pi is negative which means that the positive ones 
can take values more than unity. 

Let us discuss the condition for the system state to be separable. Accord- 
ing to Peres criterion [^Uj, the system is separable if partial transpose of the 
matrix pab (191 J) gives the positive density matrix. This condition is neces- 
sary but not sufficient. Let us discuss this condition within the framework 
of positive- map matrix representation. On the example of spin-1/2 bipatite 
system, we have shown that the map of density matrix onto its transpose can 
be included in the matrix semigroup of matrices C s . One should point out 
that this map cannot be obtained by means of the averaging with all positive 
probability distributions pk- On the other hand, it is obvious that generic 
criterion, which contains the Peres one as a partial case, can be formulated 
as follows. 

Let us map the density matrix pab of a bipartite system onto vector pab- 
Let us act on the vector pab by an arbitrary matrix, which represents the 
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positive maps in subsystems A and B. Thus we get a new vector 

pfl = (La® £b)pab. (93) 

Let us construct the density matrix p^ B using inverse map of vectors onto 
matrices. If the initial density matrix is separable, the new density matrix 
pW must be positive (and separable). 

In the case of bipartite spin- 1/2 system, by choosing La = 1 and Cb 
being matrix coinciding with the matrix g alS , we obtain the Peres criterion 
as a partial case of the criterion of separability formulated above. Thus, our 
criterion means that separable matrix keeps positivity under action of tensor 
product of two semigroups. In the case of bipartite spin- 1/2 system, the 
16x16 matrix of the semigroup tensor product is determined by 18 parame- 
ters. 

Let us discuss the positive map (J53)l which is determined by the semigroup 
for n-dimensional system. It can be realized as follows. 

The nxn Hermitian generic matrix p is mapped onto complex n 2 -vector 
p by the map described above. The complex vector p is mapped by means 
of multiplying by the unitary matrix S onto real vector p T , i.e., 

Pr = Sp, p = S- X pr. (94) 

The matrix S is composed from n unity blocks and the blocks 



where j corresponds to column and k corresponds to row in the matrix p. 
For example, in the case n = 2, one has the vector p r of the form 



Pr 



( P U \ 

V2Rep 12 
\/2Impi2 

V P22 ) 



(96) 



One has the equalities 

p1 = p e = Trp 2 . (97) 

The semigroup preserves the trace of the density matrix. Also the discrete 
transforms, which are described by the matrix with diagonal matrix blocks 
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of the form 



V 



( 1 



V o 






-1 








1 / 



(98) 



preserve positivity of the density matrix. 

For the spin case, the semigroup contains 12 parameters. 

Thus, the direct product of the semigroup (J53J) and the discrete group 
of the transform D is the positive map preserving positivity of the density 
operator. 



9 Symbols, Star-Product and Entanglement 

In this section, we describe how entangled states and separable states can 
be considered using properties of symbols and density operators of different 
kinds, e.g., from the viewpoint of Wigner function or tomogram. The general 
scheme of constructing the operator symbols is as follows [TB] . 

Given a Hilbert space H and an operator A acting on this space, let us 
suppose that we have a set of operators U (x) acting on H, a n- dimensional 
vector x = (xi,x 2 , ■ ■ ■ , x n ) labels the particular operator in the set. We 
construct the c- number function f A (x) (we call it the symbol of operator A) 
using the definition 

f A (x) = Tr [Au(x)} . (99) 

Let us suppose that relation (J9T?)) has an inverse, i.e., there exists a set of 
operators -D(x) acting on the Hilbert space such that 

A = J / i (x)£>(x) dx, Tr A = J / 1 (x)Tr£>(x) dx.. (100) 

Then, we will consider relations ()99|) and ()100|) as relations determining the 
invertible map from the operator A onto function /^(x). Multiplying both 
sides of Eq. (0) by the operator £7(x') and taking trace, one has the consis- 
tency condition satisfied for the operators f/(x') and -D(x) 

Tr [f>(x') J D(x)] = 5 (x - x) . (101) 

The consistency condition (jl01|) follows from the relation 

f A (x) = J K(x : x')f A (x')dx'. (102) 
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The kernel in (|1U2|) is equal to the standard Dirac delta-function, if the set 
of functions f A (x) is a complete set. 

In fact, we could consider relations of the form 

A -> ^(x) (103) 

and 

/ A (x) -> A. (104) 

The most important property of the map is the existence of associative prod- 
uct (star-product) of functions. 

We introduce the product (star-product) of two functions f A (x) and /^(x) 
corresponding to two operators A and B by the relationships 

Zab(x) = * /bW := Tr [iM(x)] . (105) 

Since the standard product of operators on a Hilbert space is an associative 
product, i.e., A(BC) = (AB)C, it is obvious that formula (|105|) defines an 
associative product for the functions f A (x), i.e., 

/ A (x) * (^(x) * / 6 (x)) = (^(x) * f 6 (x)) * / a (x). (106) 

Using formulas (J59j) and (jlOOj) . one can write down a composition rule 
for two symbols f A (x) and /#(x), which determines star-product of these 
symbols. The composition rule is described by the formula 

/ 4 (x) * /a(x) = / / i (x")/ i5 (x')^(x", x', x) rfx' dx". (107) 

The kernel in the integral of (|107|) is determined by the trace of product of 
the basic operators, which we use to construct the map 

K(x", x', x) = Tr [L>(x")£>(x')£(x)] . (108) 

Formula (jl08J) can be extended to the case of star-product of N symbols of 
operators Ai, A%, . . . , A^. Thus one has 



J W Ai ( Xi )W A2 ( X2 )---W An (x n ) 

xK (xi, x 2 , . . . , xat, x) cixi dx. 2 ■ ■ ■ (d$$) 
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where the kernel has the form 



if (xi,x 2 , • . • ,Xjv,x) 



Tr 



J D(x 1 )D(x 2 )--- J D(x iV )?7(x) 



110) 



Since this kernel determines the associative star-product of N symbols, it 
can be expressed in terms of the kernel of star-product of two symbols. The 
trace of an operator A N is determined by the kernel as follows: 



TyA 



N 



W A ( Xl )W A (x 2 )---W A (^ N ) 



xTr 



Z)(x 1 )/)(x 2 ) • ■ • -D(xjv) dxi dx 2 • • • dxjv- (111) 



When the operator A is a density operator of a quantum state, formula (jlllj) 
determines the generalized purity parameter of the state. When the operator 
A is equal to the product of two density operators and N = 1, formula (jlllj) 
determines the fidelity. 



10 Weyl Symbol 

In this section, we will consider a known example of the Heisenberg-Weyl- 
group representation. As operator U(x), we take the Fourier transform of 
displacement operator D(£) 

&w = /ex P - s^e) A(e»- *e, au) 

where £ is a complex number, £ = £i + z£ 2 , and the vector x = (xx, x 2 ) can 
be considered as x = (q,p), with q and p being the position and momentum. 
One can see that 

Trf>(x) = 1. 

The displacement operator may be expressed through the creation and an- 
nihilation operators in the form 

£>(£) =exp(£a t -£*a). (113) 

The displacement operator is used to create coherent states from the vacuum 
state. For the creation and annihilation operators, one has 

a= ^r- a= ^r- (114) 
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where q and p may be thought as the coordinate and momentum operators for 
the carrier space of an harmonic oscillator. The operator a and its Hermitian 
conjugate a) satisfy the boson commutation relation 

[a, a'] = 1. 

Let us introduce the Weyl symbol for an arbitrary operator A using the 
definition given by Eq. 



W A {x) = Tr 



AU(x) . (115) 



The form of operator t/(x) is given by Eq. (|112|) . One can check that Weyl 
symbols of the identity operator 1, position operator q and momentum op- 
erator p have the form 

Wi(q,p) = 1, W 4 {q,p) = q, Wp(q,p) = p. (116) 

The inverse transform, which expresses the operator A through its Weyl 
symbol, is of the form 

A= [w A (x)U(x)p. (117) 
J Ztt 

One can check that for Wj(x) = 1, formula (|117j) reproduces the identity 
operator, i.e., 

'#(x)^ = i. (118) 

Comparing ()117|) with (jlOOf) . one can see that the operator -D(x) is connected 
with £7(x) by the relationship 

fW = ^. (119) 

Let us consider now star-product of two Weyl symbols (it is usually called 
Moyal star-product). If one takes two operators A\ and A 2 , which are ex- 
pressed through Weyl symbols by formulas 

Ax = J W Al (x!)U(x>) ^ , A 2 = J iy i2 (x")t>(x") ^ , (120) 
with vectors 

x' = (x'^x'2) and x" = (x'[, x'2), 
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the operator A (product of operators A\ and A 2 ) has Weyl symbol given by 

1 

" 47T 5 



WOl(x) = Tr [AU(x)\ = J dx' dx" d 2 £d?? d 2 C 'W Ai (x')W ia (x" 

x ex P {2- 1 / 2 [(e; - + - $ + ^)(^i - z4 

+ (£1 - + *4') - (£1 + iQK - O + (6 - ^ 2 )(^i + ix 2 

-(6 + ^2)^1 -^2)] } Tr 

where £ = 6 + i£ 2 , with £' = £ + i& and £" = + 
Using properties of displacement operators 



;i2i) 



D(0^«") = ^ + Oexp(t Im(^ 



Tr 



£>(€) = ^ 2 (£), (122) 



one can get known explicit form of the kernel, which determines star-product 
of Weyl symbols. Thus we described the construction of Weyl symbols, includ- 
ing Wigner function, by means of the star-product formalism. 



11 Tomographic Representation 

In this section, we will consider an example of the probability representation 
of quantum mechanics [21] . In the probability representation of quantum me- 
chanics, the state is described by a family of probabilities [22-24]. According 
to the general scheme, one can introduce for the operator A the function 
f A (x.), where 

x = {x 1 ,x 2 ,x 3 ) = (X,/i, v), 

which we denote here as w^(X, //, v) depending on the position X and the 
parameters // and v of the reference frame 

w A (X, fx, v) = Tr [AU(x)] . (123) 

We call the function w A (X, /1, v) the tomographic symbol of the operator A. 
The operator U (x) is given by 

U(x) = U(X,»,v) = exp^(# + M)^exp^(g 2 +p 2 )^ \X)(X\ 



x ex P ^— j [Q +P ) ) exp C2P + f 

^|X)(X|C/t (124 ) 



29 



The angle 9 and parameter A in terms of the reference frame parameters are 
given by 

\i = e A cos 9, u = e~ x sin 9. 
Moreover, q and p are position and momentum operators 

q\X)=X\X) (125) 

and | X){X | is the projection density. One has the canonical transform of 
quadratures 

X = U VbV q Ul v = [iq + up, 

P = u »fr = ii£W. HEW. 

U^pU^ ^ p ^ q. 

Using the approach of j2H] one can obtain the relationship 

U (X, fi, u) = 5(X — fxq — up). 

In the case we are considering, the inverse transform determining the operator 
in terms of tomogram [see Eq. (|1(J(J|) ] will be of the form 

A = J w A (X, fjt, v)D(X, fi, v) dX dyL du, (126) 

where ^ 

-D(x) = D(X, fM, v) — — exp (iX — iup — i/xq) , (127) 

i.e., 

D{X,n,v) = -±-exp(iX)D(t(v ll v)). (128) 
The unitary displacement operator in (J128j) reads now 

= ex P (t(p, ~ £*G", v)a) , 

where £(//, u) = £1+^2, w ^ n Ci = R- e (£) = v/\/2 and £ 2 = Im (£) = — /V\/2. 

Trace of the above operator which provides the kernel determining the 
trace of an arbitrary operator in the tomographic representation reads 

Tr/)(x) = e lX 5{fi)5{u). 
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The function w A (X, //, v) satisfies the relation 



w A (AX, A//, Az/) = T—r ^(X, fj,, v). 



(129) 



This means that the tomographic symbols of operators are homogeneous 
functions of three variables. 

If one takes two operators A\ and A 2 , which are expressed through the 
corresponding functions by the formulas 



A l = j w Ai (X' >' ',u')D(X' >' \v')dX' 'dfj,' 'dv' ', 
A 2 = J w A2 (X", //, u")D(X", //, v")dX" dfi" dv", 



(130) 



and A denotes the product of A\ and A 2 , then the function w A (X, /i, v), which 
corresponds to A, is star-product of functions w A (X, /i, v) and w A (X, /i, v). 
Thus this product 

w A (X, /i, v) = w Ai (X, fi, v) * w A , 2 (X, /i, u) 

reads 

w A {X,n,v) = J Wii (x") Wi2 (x / )X(x",x',x)rfx"rfx', (131) 
with kernel given by 

K(x", x', x) = Tr [£)(X", //, u")D(X', //, i/)U(X, fi, v)] . (132) 
The explicit form of the kernel reads 

K(X 1 ,fx 1 , v u X 2 , /i 2 , ^2, Xp, v) 



4tt 2 

1 - VI - 



The kernel for star-product of N operators is 

K (X 1 ,/j, 1 , v u X 2 , nz, v 2 ,..., X N , fj, N , i/jv, X, n, v) 

N 



exp ( -{ (z/!/i 2 - U2H1) + 2Xi + 2X 2 

xlV.133) 



(1/1 + i/ 2 ) H (A*i + A*2j 



5(/iEf =1 ^-^Ef =1 ^) /i J " . . " v 

ex P I o \ 2^ {l/L t'j ~ ^/"■ • ) • -2^ X j 



(2 

i-Vr-W^. | + 



fe<j=i 



xl).(134) 
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The above kernel can be expressed in terms of the kernel determining star- 
product of two operators. 



12 Multipartite Systems 

Let us assume that for multimode (iV-mode) system one has 

U(y) = Uu(^ k) ), (135) 
k=i 

N 

D(y) = J] D (fW) , (136) 



fc=i 
where 

i7=(4 1) ,4 1) ,...,xS ) ,x?\4 2) ,...,#). (137) 

This means that symbol of density operator of the composite system reads 

m = Tr[pI[U(x {k) )}- (138) 
k=i 

The inverse transform reads 

. N N m 

P = / dyf p (y) I] D(x^), dy=UU dxfl (139) 

k=l k=ls=l 

If the symbol corresponds to a system Wigner function W(q,p), the operator 
U(x^), where = (q k ,Pk), was discussed above. It has the form 

U(x^) = V(a k )(-l) a l ak V(-a k ), (140) 

where 

ak = 71 ( qk+ipk ^ V ^ = e afc °*~ a *° fc , [a k , a m ] = 0, [a k , 4j = 5 km . 

(141) 

The operator 

D(x k ) = -U(x {k) ). (142) 

The tomographic symbols are also defined by analogous formulas with specific 
x and operators U(x k ), D{x k ). 
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Now we formulate properties of symbols in the case of entangled and 
separable states, respectively. 

Given a composite m-partite system with density operator p. 

If the nonnegative operator can be presented in the form of 'probabilistic 
sum' 

P = 5>O0p?° ® P? 2) ® • ■ • ® P? m \ (143) 

with positive probability distribution function V(z), where components of z 
can be either discrete or continuous, we call the state 'separable state'. This 
means that the symbol of the state can be presented in the form 

m 

f P (y)=T,n?)Ilrt ak) (xk^- ( 144 ) 

z k=l 

For example, the Wigner function of separable state of bipartite system has 
the form 

W( qi ,q 2 , Pl ,p 2 ) =EP(^ (1) (?i.Pi.^ (2) (92.P2,4 (145) 

z 

Analogous formula can be written for the tomogram of separable state. 

13 Spin Tomography 

Below we concentrate on bipartite spin systems. 

The tomographic probability (spin tomogram) completely determines the 
density matrix of a spin state. It has been introduced in [TD| ITU ITS]. 

The tomographic probability for spin-j state is defined via the density 
matrix by the formula 

(jm | D\g)pD{g) \ jm) = W®(m, 0), m = -j, -j + 1, . . . , j, (146) 

where D(g) is the matrix of S'?7(2)-group representation depending on the 
group element g determined by three Euler angles. The set of the tomogram 
values for each is an overcomplete set. We need only finite number of 
independent locations which will give information on the density matrix of 
the spin state. Due to structure of the formula, there are only two Euler 
angles involved. They are combined into the unit vector 

= (cos 4> sin d, sin sin i?, cost?). (147) 
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This is the map from S 3 to S 2 . 

The physical meaning of the probability W(m, 0) is the following. 

It is the probability to find, in the state with the density matrix p, the spin 
projection on direction equal to m. For bipartite system, the tomogram is 
defined as follows: 

W(mim 2 1 2 ) = (jirnij 2 m 2 \ D ] \gi)D ] \g 2 )pD{g x )D{g 2 ) \ j\m\j 2 m 2 ). 

(148) 

It completely determines the density matrix p. It has the meaning of joint 
probability distribution for spin j\ and j 2 projections mi and m 2 on directions 
0i and O2. Since the map p # W is linear and invertible, the definition of 
separable system can be rewritten in the following form of decomposition of 
the joint probability into sum of products (of factorized probabilities): 

Wim^oA) =J2PkW ik) (m 1 6 1 )W ik) (m 2 6 2 ). (149) 
k 

This form can be considered to formulate the criterion of separability of the 
two spin state. 

The state is separable iff the tomogram can be written in the form ()149|) 
with J2kPk = 1, Pk > 0. It seems that we simply use the definition but, 
in fact, we cast the problem of separability into the form of property of 
the positive joint probability distribution of two random variables. This 
is area of probability theory and one can use results and theorems on the 
joint probability distributions. If one does not use any theorem, one has to 
study solvability of relation (|149|) considered as the equation for unknown 
probability distribution p& and unknown probability functions W^ k \mJ)i) 
and W^ k \m 2 6 2 ). 

14 Example of Spin- 1/2 Bipartite System 

For spin-1/2 state, the generic density matrix can be presented in the form 

p=^(l + a-n), n = (n 1 ,n 2 ,n 3 ), (150) 

where a are Pauli matrices and n 2 < 1, with vector n for a pure state being 
unit vector. This decomposition means that we use as basis in 4-dimensional 
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vector space the vectors corresponding to Pauli matrices and unit matrix, 
i.e., 
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V o J 
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-1 J 




V i / 



(151) 



The density matrix vector 

/ Pn \ 

Pl2 
P21 
V P22 J 

is decomposed in terms of the basis vectors 



P 



(152) 



?= 5 (i + 



+ ^3^) • 



(153) 



It means that tomogram of spin-1/2 state can be given in the form 



Wi-,0 



1 n- 

2 + — 



1 ra-0 

2 2~~ 



(154) 



Inserting these probability values into relation 
get the relationships: 



for each value of k we 



1 1 

2'2' 



0i,0 2 ) = \+\ {y,P k n^A+\ (j2p k n*Xo 2 +J2Pk (n k • $i) (n 8 k ■ 0, 

(155) 

,0i,0 2 ) = ^ (Ep^)-Oi-^ (£^)-02-£p* (n fc ■ Ox) (n£ ■ 2 ) 

(156) 



1 1 

2'~2 : 



One has the normalization property 



1/2 



E ^(mimsOiOs) = 1. 

mi, 7712 = — 1/2 



(158) 
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One easily gets 



w Q, i o 1; o 2 ) + w (i -± o 1; o 2 ) = i + (l>"*) • Oi- 

This means that derivative in Oi on the left-hand side gives 



_d_ 

Analogously 

_d_ 

W 2 



W 



1 u . 



1 1 



w 



T 2 



w o.o,0i,02 )+w [--,-, 0x,0 



1 1 



n 



(*) 



(159) 



(160) 



(161) 



Taking the sum of (J156)) and (j!57)l ) one sees that 



1 d d 



W 



• o L . Go ) -1 ir(-H,Oi,o 2 



-I^PfcK) 

k 



(162) 



Since we look for solution where p k > 0, we can introduce 

N k = VP~A, N^ = ^p- k n^. (163) 
This means that the derivative in ()162|) can be presented as tensor 

= £(^ < (tf£ r) ) i . (164) 



One has 

(165) 
(166) 

k k 

The conditions of solvability of the obtained equations is a criterion for sep- 
arability or entanglement of bipartite quantum spin state. 
For Werner state j2E] with the density matrix 



k k 



Pab 



f (l+p)/A p/2 \ 

(l-p)/4 

(l-p)/4 

V p/2 (l+p)/Aj 



1/10 



(167) 
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one can reconstruct known results that for p < 1/3 the state is separable 
and for p > 1/3 the state is entangled, since in the decomposition of density 
operator in the form ()149|) the state 



Po = 7 



( 1 \ 

10 

10 

V o o o i j 



(168) 



has the weight po = (1 — 3p)/4. 

For p > 1/3, the coefficient p becomes negative. 

There is some extension of the presented consideration. 

Let us consider the state with the density matrix (nonnegative and Her- 
mitian) 



/ Rn 



V R21 





Pn 
P21 





P12 
P22 







R\2 




R22 



\ 



Tr p = 1 . 



(169) 



Using procedure of mapping the matrix onto vector p and applying to the 
vector nonlocal linear transform corresponding to Peres partial transpose and 
making inverse map of the transformed vector onto the matrix, we obtain 



( Rn 

Pn 

R21 

V P21 



P12 \ 

R12 

P22 

R22 J 



(170) 



In the case of separable matrix p, the matrix p m is nonnegative matrix. 
Calculating eigenvalues of p m and applying condition of their positivity, we 
get 



R11R22 > \Pl2\ 



P11P22 > \R 



12 



(171) 



Violation of these inequalities gives a signal that p is entangled. For Werner 
state (ITKTf) . Eq. (fTTTl) means 



1 +p > 0, 



1 — p > 2p, 



(172) 



which recovers the condition of separability p < 1/3 mentioned above. 

The joint probability distribution ()148J) of separable state is positive af- 
ter making the local and nonlocal (Peres-like) transforms connected with 
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positive map semigroup. But for entangled state, function (|148|) can take 
negative values after making this map in the function. This is a criterion of 
entanglement in terms of tomogram of the state of multiparticle system. 



15 Dynamical Map and Purification 

In this section, we consider connection of positive maps with purification 
procedure. In fact, formula 



P ft = YsPkUkpUl, 



(173) 



where Uk are unitary operators, can be considered in the form 

P^p' = J2Pkpk, Pk>0, Y.Pk = 1 - ( 174 ) 



Here the density operators pu read 

Pk = U k pW k . 



(175) 



This form is the form of probabilistic addition. This mixture of density 
operators can be purified 



p'^p" = N 



'TipkPopjPo 



;i76) 



where P is a fiducial projector and 



pkPopj 



Tr p k P pjP 



(177) 



The map ()173j) could be interpreted as the evolution in time of the initial 
matrix p Q considering unitary operators Ukif) depending on time. Thus one 
has 

Po - p(t) = J2PkU k (t) P oUt(t). (178) 

k 

In this case, the purification procedure provides the dynamical map of a pure 

state 

\ipa)(ipo\^\m)(*P(t)\, (179) 
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where | ip(t)) obeys to a nonlinear equation and, in the general case, this 
equation is not differential equation in time variable like the Schrodinger 
equation. 

For a specific case, the evolution (j!78J) can be described by semigroup. 
The density matrix ()178|) obeys to first-order differential equation in time for 
this case [27-29]. 

The reason why there is no differential equation in time for generic case 
is due to the absence of the property 

PiAh) =Y, K T^)P™{tx), (180) 

mn 

where the kernel of evolution operator satisfies 

K% n (t 3 ,t 2 )KZ(t2,h) = Kff(t 3 ,h). (181) 

It means that trajectory (curve) is not determined by differential equation 
in time. 

Thus via purification procedure and dynamical map of the density matrix 
we get the dynamical map of a pure state (nonlinear dynamical map). This 
map can be used in nonlinear models of quantum motion. 



16 Density Matrix and Real Quadratic Forms 

It is convenient to associate the Hermitian nonnegative nxn density ma- 
trix p with real quadratic form determined by the real matrix D using the 
relationships 

r = P Af, iR = P ~^f> rt = r > R l = ~R, (182) 

and 

The quadratic form is the scalar function (homogeneous polynomial of second 
order) 

f(x,y,r,R) = (x,y)D^^ > 0. (184) 

Nonnegativity of f(x, y, r, Flj takes place for all nonnegative Hermitian ma- 
trices p. 
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In fact, for complex vectors, 

z = x + iy, (185) 

one has 



z*pz = f(x,y,r,R). (186) 
All real transforms of the form 

D -> D' = ADA T (187) 

keep the quadratic form nonnegative. 

Let us consider the real 2nx2n matrix A given in block form 

A -{'dy (i88 » 

where a, b, c, and d are real nxn matrices. 
If the matrix ()187j) has the form 

D ' = ( R* / ) ' < 189 > 

one has the relationships: 

ara* + fciZV + aRtf + brtf = r', (190) 

crc* + dRV + cRd* + drd* = r'\ (191) 

arc* + oi?V + aRd t + brd* = R' , (192) 

era* + dRV + ciZo* + drtf = R'\ (193) 

Also R n = -R', r !t = r'. 
For the case 

b = c = 0, (194) 

one has possible solutions 

a = d (195) 

and 

a = -d. (196) 
Solution (|196|) for a = 1 describes a Peres-like transposition of the matrix p. 
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For the case a = d = 0, one has possible solutions 



b = c 



(197) 



and 



b = —c. 



(198) 



In the case ()195|) . one has 



r' = ara*, R' = a-Ra*, 
r' + iR' = a(r + iR)a\ 



(199) 
(200) 



In the case (|196|) . one has 



r' = ara 1 , R' = —aRa 1 , 
r' + iR' = a(r -iR)a\ 



(201) 
(202) 



Thus for block diagonal matrices A the possible transforms of the initial 
density matrix matrix p have the form 



In the vector form p <-» p, the transform ()203|) is described by superoperators 



Obviously one can apply the averaging procedure to get the matrix (a <g> 
a). It is a partial case of the general transform (v <8> v) for real v. Here 
superpoperator 1} makes from the vector p ^ p the vector p t <-» p l , where 
p t is transposed density matrix p. 

The superoperator V" in the case n = 2 is described by the matrix g a @. 
The solutions (|197|) and ()198|) provide analogous transforms ()203|) and (j204J) 
with replacement a —> b. 

For n = 2, the choice a = o"i gives 




(203) 




(204) 



(205) 



which is exactly Peres transpose transform. 
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For n — 3, the choice 




(206) 



provides two transforms of the Hermitian density matrix 



and 



( Pll Pl2 Pl3 
P21 P22 P23 
\ P31 P32 P33 



P -»• P 



( Pll Pl2 
P21 P22 
V -P31 -P32 



-Pl3 N 
-P23 
P33 / 



(207) 




(208) 



Obviously, the unitary transform u of the form 

p — > upv) 

does not change the nonnegative eigenvalues of the density operator but this 
transform differs from the transforms discussed above, e.g., if |det a\ ^ 1, the 
above transforms do not preserve the determinant of the density matrix. 

The given construction can provide also noncompletely positive map. For 
example, if the transform of Hermitian nonnegative 3x3 matrix pjj, is de- 
scribed by the formula 

Pjk —> Pjfe( cosh 2 9 cos(9j — 6 k ) — sinh 2 65jkj 

(there is no sum over j, k) , it corresponds to applying to vector p the matrix 
(a ® a) with real matrices a. The formula can be used also for arbitrary 
integer n. But averaging is done using quasidistribution (not the probability 
distribution). This means that in sum Sfc £ fc a fc® a fc the numbers e k take both 
positive and negative values ±1. It is another example of noncompletely 
positive map. The example of Peres transpose transform discussed in the 
previous sections for the case n = 2 belongs also to the case of positive but 
not completely positive maps. 
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17 Tomogram of the Group U(ri) 

In order to formulate a criterion of separability for a bipartite spin system 
with spin ji and j 2 , we introduce the tomogram w(l,m, g( n >) for the group 
U(n), where 

n = n x n 2 , n x = 2j 1 + 1, n 2 = 2j 2 + 1, 

and are parameters of the group element. Vectors I and in label a basis 
I Z,m) of the fundamental representation of the group U(n). For example, 
since this representation is irreducible, being reduced to representation of 
U{n\) <8> t/(^2)-subgroup of the group U(n), the basis can be chosen as the 
product of basis vectors: 

\ji,mi) \h,m 2 ) ^\ji,j 2 ,mi,m 2 ). (209) 

Due to irreducibility of this representation of the group U{n) and its sub- 
group, there exists a unitary transform u l ^ 2mim2 \ l,m) such that 

I ji,j2,m h m 2 ) =J2 u hhmim 2 \l,fn), (210) 

im 

I I™) = (""^iSamjma I 3u h, m u m 2 ). (211) 

m\rn,2 

One can define the U(n) -tomogram for a Hermitian nonnegative nxn density 
matrix p, which belongs to Lie algebra of the group U(n), by a generic 
formula: 

w(tm,g (n) ) = (T,m I U\g {n) )pU(g {n) ) \ l,m). (212) 

Formula ()212|) defines the tomogram in basis | l,rh). 

Now let us define the U(n) -tomogram using basis | ji, j 2 , mi, m 2 ), i.e., 

w { j^\m 1 ,m 2 ,g {n) ) = (ji,j2,m 1 ,m 2 \ U\g {n) )pU(g {n) ) \ ji,j 2 ,m x ,m 2 ). 

(213) 

This tomogram is spin-tomogram P31 for G U(2)®U(2) sub group of the 
group U(n). Properties of this tomogram follow from its meaning to be joint 
probability distribution of two random spin projections mi, m 2 depending on 
g^ parameters. 

One has normalization condition 

w^ h \mi,m 2 ,g^) = 1. (214) 

mi ,m2 
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Also all the probabilities are nonnegative, i.e., 

w Ul ' j2 \mi,m 2 ,g {n) ) > 0. (215) 



Due to this, one has 

For spin-tomogram, 
and 



\w in ' j2 \m h m 2 ,g {n) )\ = l. (216) 



mi, m,2 



9 {n) ^(d 1 ,d 2 ) (217) 



w (jl ' j2 \mi, m 2 , g {n) ) -> w(m 1 ,m 2 , O u 2 ). (218) 

The separability and entanglement condition discussed in the previous 
section for bipartite spin-tomogram can be considered also from the viewpoint 
of the properties of £/(n)-tomogram. If the two-spin nxn density matrix p 
is separable, it keeps to be separable under action of generic positive map of 
the subsystem density matrices. This map can be described as follows. 

Let p to be mapped onto vector p with n 2 components. The components 
are simply ordered rows of the matrix p, i.e., 

P = (pll ,Pl2, ••• , Pin, P21, P22, ••• ,Pnn,)- (219) 

The n 2 xn 2 matrix L is taken in the form 

L = Y.PsL^®L^\ p s >0, $> s = l, (220) 

s s 

where n\Xni matrix and n 2 xn 2 matrix describe the positive maps 
of density matrices of spin-jx and spin-j 2 subsystems, respectively. We map 
vector p onto vector pi 

PL = Lp (221) 

and construct the n 2 xn 2 matrix p L , which corresponds to the vector p L . 
Then we consider U(n) -tomogram of the matrix pi, i.e., 

w ( l un) {m u m 2 ,g {n) ) = (ji,j 2 ,m 1 ,m 2 | U\g {n) p L U(g {n) ) \ j u j 2 ,m h m 2 }. 

(222) 

Using this tomogram we introduce the function 



F(gW,L)= E \w { r 2 \mum 2 ,g^)\. (223) 
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For separable states, this function does not depend on the £/(n)-group pa- 
rameter and positive-map matrix elements of the matrix L. 

For normalized density matrix p of the bipartite spin-system, this function 
reads 

FfoW L) = l. (224) 

For entangled states, this function depends on g( n ' and L and it is not equal 
to unity. This property can be chosen as necessary and sufficient condition 
for separability of bipartite spin-states. In fact, the formulated approach can 
be extended to multipartite systems too. The generalization is as follows. 

Given N spin-systems with spins ji, j 2 , ■ ■ ■ ,Jn- Let us consider the group 
U(n) with 

N 

n = n n k , n k = 2j k + 1. (225) 
k=i 

Let us introduce basis 

N 

I m) = II I jkmk} (226) 

k=l 

in the linear space of the fundamental representation of the group U(n). We 
define now [/(n)-tomogram of a state with n 2 xn 2 matrix p: 

w p (m,g {n) ) = (m | U\g {n) )pU{g^) \ m). (227) 

For positive Hermitian matrix p with Trp = 1, we formulate a criterion of 
separability as follows. 

Let the map matrix L to be of the form 

N 

L = E^(ll®^ fc) )' P»>0> £Ps = l, (228) 

s k=l s 

where is positive-map matrix of the density matrix of fcth spin subsystem. 
We construct the matrix pi as in the case of bipartite system using the matrix 
L. The function 

F(^,L) = ]TKJm,^))|>l (229) 

m 

is equal to unity for separable state and it depends on the matrix L and U{n)- 
parameters g( n ' for entangled states. This criterion can be applied also in the 
case of continuous variables, e.g., for Gaussian states of photons. Function 
(J229|) can provide the measure of entanglement. Thus one can use maximum 
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value (or a mean value) of this function as a characteristic of entanglement. 
In fact, the separability criterion is related to the following positivity criterion 
of finite or infinite matrix A. The matrix A is positive iff the sum of moduli 
of diagonal matrix elements of the matrix UAW is equal to positive trace of 
the matrix A for arbitrary unitary matrix U. 

18 Conclusions 

To conclude, we formulated the notion of separability and entanglement as 
a criterion for joint tomographic probability of subsystem states to be repre- 
sented in the specific form of sum of products of tomograms of the subsys- 
tems. 

We have shown that the positive map of density matrix of multiparticle 
system expressed in terms of superoperator acting in Lie algebra (adjoint 
representation) of unitary group can be considered as a semigroup, which 
contains all local unitary transforms acting in subspaces corresponding to 
the subsystem states. 

The set of separable states is shown to be invariant under action of this 
group. 

The intrinsic measure of entanglement is shown to be invariant under 
action of the local group. 

The formalism of vectors representing the matrices is convenient tool for 
the consideration. We introduced unitary spin tomogram and formulated 
necessary and sufficient condition of entanglement. 
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